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Abstract 

The equations for fully compressible rotating magnetoconvection are numerically 
solved in a Cartesian box assuming conditions roughly suitable for the geodynamo. 
The mean electromotive force describing the generation of mean magnetic flux by 
convective turbulence in the rotating fluid is directly calculated from the simula- 
tions, and the corresponding a-coemcients are derived. Due to the very weak density 
stratification the a-effect changes its sign in the middle of the box. It is positive at 
the top and negative at the bottom of the convection zone. For strong magnetic fields 
we also find a clear downward advection of the mean magnetic field. Both of the 
simulated effects have been predicted by quasi-linear computations (Soward, 1979; 
Kitchatinov and Riidiger, 1992). Finally, the possible connection of the obtained 
profiles of the EMF with mean-field models of oscillating a 2 -dynamos is discussed. 
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1 Introduction 



Convective motions in the fluid outer core influenced by rotation and magnetic 
fields are able to maintain the Earth's magnetic field for long times interrupted 
by occasionally occurring reversals of the dominating dipole component. The 
equations describing the physical processes in the fluid outer part of the 
Earth 's core are rather stiff f Br aginskv and Robert], Il99i I iRoberts and Glatzmaierl . 
2000) showing a broad range of timescales on which the characteristic behavior 
is observable. The resulting and dominating advective timescale r adv = d/v! 
leads to a very short coherence length of the physical quantities like veloc- 
ity or magnetic field which, as a consequence, requires a very high spatial 
resolutio n in numerical sim ulations to include the effects of small-scale tur- 
bulence (|Hollerbachl . 120031b In case of a conducting and rotating fluid with 
large-scale density stratification a convectively driven turbulence can generate 



Preprint submitted to Elsevier Science 



2 February 2008 



a mean electromotive force (EMF) parallel to the m ean magnetic field - a 
process known as a-effect f Krause and Radler . 1980). But the Earth's fluid 



outer core, though a rapid rotator, is rather weakly stratified so that it re- 
mains unclear whether the amplitude of the a-effect is sufficient to maintain 
dynamo action. 

Present-date global numerical simulations are unable to include small-scale 
motions because of computational restrictions. Nevertheless, they show quite 
satisfying results in the sense th at they are able to reproduce several features of 
the observed geomagnetic field (iGlatzmaier and RobertsL 19961: Christensen et a! 



1999; iKuang and Bloxhaml . ll999l ). In virtue of these results it seems to be justi- 



fied to ignore small-scale fluctuations but, on the other hand, there still remain 
many open questions when comparing such sim ulations with their restricted 



param eter range with observational data (see e.g. lDormv. Valet, and Courtillot . 
200fJ. 



Background of our calculations is the idea to describe a dynamo by an induc- 
tion equation for the mean magnetic field (B) with a prescribed productive 
term due to the small-scale turbulence that ensures the existence of dynamo 
action. 

In mean-field dynamo theory the magnetic field and the velocity are split up 
in a mean part, (B) and (u), and a fluctuating component, B' and u' . The 
time behavior of the mean magnetic field (B) = B — B' is described by 

d t (B) = Vx ((it) x (B) + £—r]V x (B)) (1) 

with the mean electromotive force S = (u' x B') which is usually expressed 
by 

Si = («' x B') i = aij (Bj) + (3 ijk d k (B 3 ) . (2) 



The tensor correlates the turbulent EMF due to small-scale motions with 
the large-scale magnetic field including the effects of anisotropy. The tensor 
flijk is related to the turbulent diffusivity t]t by (3ijk = rj^eijk- The easiest way 
to look for dynamo action is to solve equation (1) with zero mean flow (u) and 
with an EMF taken from equation (2) without turbulent diffusivity (/%& = 0). 

Oscillating solutions f or so- called a 2 -dynamos have recently been presented by 
Stefani and Gerbethl ((2003) . They determined a radial profile of a spherically 



symmetric and isotropic a under the conditions that the dynamo mode with 
the lowest eigenvalue is an oscillating solution and all stationary modes (that 
usually dominate a 2 -models) are damped. This constraint leads to an a-effect 
with characteristic zeros in the radial pr ofile. Calculations w ith a uniform 
radial a-coefficient have been performed by Riidiger et al. ( 2003| ) who, instead, 
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included a latitudinal variation of a and effects of anisotropy. They found 
oscillating a 2 -dynamos only for exotic exceptions. 



Hovng et al.l ([20011 ) analyzed the behavior of a geo dynamo- like aJ?-dynamo 
model which produces an axial dipole with reversals that are induced by a 
fluctuating a. The amplitude of the fundamental dipole mode behaves as a 
damped particle under the influence of a random force in a bistable potential. 
Here, in contrast to that, a configuration without any differential rotation and, 
hence, without i2-effect is considered. 

The purpose of the present paper is to calculate the turbulent EMF directly 
from numerical solutions of the full set of nonlinear MHD-equations for a con- 
vectively driven turbulent flow under the influence of rotation and subject to 
an imposed magnetic field. Our aim is to get a representation for the EMF 
based on simulations under conditions characteristic for the geodynamo. The 
results principally include effects of anisotropy of the a-effect, its radial de- 
pendence and its quenching properties. The derived a-coefficients will serve 
as some input data for future mean-field a 2 -dynamo calculations. A simplified 
geometry of a Cartesian box representing a small part of a rotating spherical 
shell provides the ability to examine the small-scale behavior of the fluid mo- 
tions and magnetic field and allows to consider effects of the turbulence that 
are neglected in global simulations. 



2 The model 



2. 1 General properties 



Our m odel i s an ada p tion o f configurations used, for example, bv lOssendriiver et al 



f)20mh and IZiederl (|2002h who examined rotating magnetoconvection in a 



box suitable for the solar convection zone. Figure 1 shows a sketch of the 
box placed somewhere on a spherical shell at some latitude 9. The coordi- 
nate system is chosen such that the unit vectors x,y, z form a right-handed 
corotating system with x pointing towards the equator, y pointing in the 
toroidal direction (from west to east) and z pointing from the bottom to the 
top of the box. Translating this Cartesian system into global spherical coor- 
dinates, z represents the radial direction f directed from inside to outside, 
y the azimuthal direction <p and x the meridional direction 6, respectively. 
The angular velocity ft in the local box coordinate system is then given by 
£2 = —Q smOx + Q cosOz where is the angular velocity of the rotating 
spherical shell. 

We try to construct a simple model that, at least, roughly represents the 
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Fig. 1. Model box being part of a rotating spherical shell at latitude 6. 

conditions in the Earth's fluid interior. Our Cartesian box model consists of 
one convectively instable layer with a weak density stratification. The param- 
eters have been chosen in a way that, ultimately, a low Mach number flow 
with Ma ~ O(10~ 2 ) results and compressibility effects, though existing, are 
rather small. We further restrict our computations to a rapidly rotating box 
expressed by a Rossby number Ro = v! /2fld <C 1. Our principal interest 
focuses on the presence of strong magnetic fields with significant dynamical 
influence on the flow. To investigate the transition from the weak field case to 
the strong field case the strength of the imposed magnetic field is successively 
increased covering a wide range of magnitudes. Here, we present results ob- 
tained for a magnetic field applied in y-direction corresponding to a toroidal 
field in spherical coordinates i.e. the production of a poloidal field from a 
toroidal field via the a-coefficient a yy is examined. Imposing other field com- 
ponents will be the subject of subsequent studies. In the strong field case this 
configuration ensures that the inertial and viscous forces are negligible and 
the main balance between the forces governing the magnetoconvection state 
is given by the Corioli s force and the Lorentz force as it is supposed to be the 
case in the fluid core (|Hollerbachl . Il99fit IZhang and Schubert! lioooh . 



2.2 Equations 



The MHD-equations for a rotating fluid including the effects of thermal con- 
duction, compressibility, viscous friction and losses due to magnet i c diffu sivity 
are solved numerically using the code NIRVANA (Ziegler, 19981 19991 ). The 
equations in the local corotating system are 



d t p=-V ■ (pu) 



(3) 



d t (pu) = -V- (puu)-VP+V-a+pg-2pf2xu+—(VxB)xB (4) 

Ho 
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d t e = -V ■ (eu) - PX7-u + ao\7u + — \VxB\ 2 + V-(xVT) (5) 
a t B = Vx(uxB-? ? VxB) (6) 

with the density p, velocity u, pressure P, magnetic flux density B, tempera- 
ture T and the thermal energy density e. We assume a constant gravitational 
field g = —gz within the domain. The viscous stress tensor a is given by 
<y%3 = vp(dui/dxj + duj/dxi — 2 / , zV-u5ij). v denotes the kinematic viscosity 
and x the thermal conductivity coefficient. The values of x, the dynamic vis- 
cosity z/dyn — V P and the magnetic diffusivity rj are constant over the box vol- 
ume. An ideal gas equation of state is assumed with P = (7— l)e = k(mp)~ l pT 
where k is the Boltzmann constant, m the atomic mass unit, p the mean 
molecular weight (p = 1 for all runs) and 7 = C P /C V = 5/3 is the ra- 
tio of the specific heats. The permeability p is given by the vacuum value 
p = 4tt x lO^VsA^m- 1 . 



2.3 The initial state 



From the equation of state and the condition for hydrostatic equilibrium, 
d z P = —pg, together with the assumption of a polytropic temperature dis- 
tribution, T = T (p/p ) r , the initial density distribution can be calculated 



as 



l/r 



p(z)=p U + ^(d-z) (7) 



where d stands for the vertical box extension and the polytropic index V is 
given by T = In (1 + dd z T/T ) / ln£. The stratification index £ = Pbot/Pto P , the 
temperature T and the global temperature gradient d z T are prescribed input 
parameters whose values are given below. The subscript refers to values 
taken at the top boundary of the box. 

The gravitational acceleration can be calculated from the hydrostatic equilib- 
rium condition and the initial density distribution (7) and is given by 

9 = ^—dzT. (8) 
T mp 



To obtain a convectively unstable state the condition T > 7 — 1 must be 
fulfilled. In fact, T = 7 — 1 leads to a Rayleigh number 

XV \ CpJ 
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with Cp = k(mn) x 7(7 — 1) 1 the specific heat at constant pressure. Here, 
parameters are chosen such that T > 7 — 1. 



2.4 Boundary Conditions 



All quantities are subject to periodic boundary conditions in the horizontal di- 
rections. At the top and at the bottom of the computational domain constant 
values for density and temperature are imposed. The vertical boundary con- 
dition for the magnetic field is a perfect conductor condition, and a stress-free 
boundary condition is adopted for the horizontal components of the velocity 
u x and u y . Impermeable box walls at the top and at the bottom lead to a 
vanishing u z at the vertical boundaries. Table 1 summarizes these conditions 
and gives the initial values for density and temperature which describe the 
overall stratification and the global temperature gradient. 



Table 1 
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2. 5 Input parameters 



All input quantities are measured at the top of the box. This makes sense 
since the density variation with depth is negligible and the temperature varies 
only by a factor of 2. The parameters is, x, V and J? are calculated from the 
Rayleigh number Ra as defined above, the Prandtl number Pr = vpCp/x, the 
magnetic Prandtl number Pm — vjr\ and the Taylor number Ta = 4J? 2 rf 4 /z/ 2 . 
The basic parameter set used for all simulations that are presented in this 
paper is given by Ra = 10 6 , Pr = 0.5, Pm = 0.5 and Ta = 10 7 . 

The Elsasser number 

A = ^F (10) 
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serves as an input parameter for the magnitude of the imposed magnetic field 
whose influence is investigated by varying A from 10~ 2 to 10 3 covering the 
full range from weak fields to very strong fields. The box with an aspect ratio 
8:8:1 is placed at a latitudinal angle of 45° on the northern hemisphere of the 
rotating spherical shell and a standard resolution of 100 x 100 x 80 grid points 
is used in all calculations. For all simulations temperature and density at the 
top of the box are scaled to unity, as it is the case for the global temperature 
gradient d z T and the box height d. A stratification index of £ = 1.1 is used. 



3 Results 

3.1 General properties and energetics 

At first, a non- rotating non- magnetic convection model is computed and the 
resulting statistically steady state is used as initial condition for the full prob- 
lem of rotating magnetoconvection. Typical values for the turbulent velocity 
u' and the turbulent magnetic field B' are obtained by an averaging proce- 
dure that includes the whole box volume. In the following, volume averages 
are indicated by double brackets, ((•)), whereas horizontal averages are de- 
noted by single brackets, (•). As root mean square value of fluctuations we 
define ((f 2 » = (N x N y N z )- 1 £ (f ijk - {f) k f. Here, N x (N y , N z ) denotes the 

number of grid cells in x, (y, z) direction and fyk — (f) k is the deviation of 
the fluctuating quantity at a certain grid cell labeled ijk from its horizontal 
average. Note that due to the horizontal averaging procedure and the periodic 
horizontal boundary conditions the mean quantities have no dependence on x 
or y. 

Time averages are labeled by / and are computed only over time intervals 
that show no significant change in the average itself. In the following, time 
is measured relative to the turnover time given by r adv = d/u rms , where 

u Tms = <JJJu^yj. All time averages are calculated over a time range of at 

least 20r a( jv starting at a certain time after the effects of magnetic field and 
rotation have been introduced and ((u' 2 ^ has reached a new statistically 
steady state. A comparison with long-term computations shows that this is a 
sufficient timespan in order to obtain meaningful results. The longest run has 
been performed for A = 4 - namely more than 110 turnover times - which 
corresponds to about two magnetic diffusion times r v = d 2 /rj. Note that r a d v 
evolves as a part of the solution and depends on the imposed magnetic field 
whereas t v can be obtained from the input parameters. 

Figure 2 shows the temporal behavior of the kinetic energies E£?' z = J dVpu 2 x y /2 
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Fig. 2. Time dependence of the kinetic energy and the magnetic energy, 
x-component (solid), y-component (dotted), z-component (dashed). A = 1. 

(left) and the magnetic energies E^gg = J dV(2fj, )~ 1 Bl y z (right) for the sim- 
ulation with A = 1. The first part of the simulation up to t ~ 17r a H v corre- 
sponds to thermal convection without rotation and without magnetic field. 
At t ~ 17r a dv the effects of rotation and magnetic field are added seen by 
a significant drop in the kinetic energies (left panel of figure 2) and a sharp 
increase of magnetic energies (right panel of figure 2). After a short tran- 
sition phase the energy of the different components remains approximately 
constant for the remainder of time. A large-scale magnetic field in x-direction 
establishes during the timespan t = 17r a dv-20r a d v (see also Figure 5 below) 
which is associated with a remarkable amount of magnetic energy stored in 
the ^-component (solid line in the right panel of figure 2). 

Since the induction equation (1) for the mean magnetic field gives d t (B z ) = 
no (B z ) can evolve during the simulations. Therefore, the magnetic energy in 
the vertical field component results from the fluctuating component B' z . The 
time dependence of the energies is qualitatively similar for all runs. Major 
differences appear in the amount of quenching after the effects of rotation and 
magnetic field have abruptly been introduced. 



Turbulent quontities Energy 




Fig. 3. Left: Quenching of turbulent velocity and (normalized) magnetic field fluc- 
tuations. Right: Ratio of the magnetic energy to the kinetic energy as a function of 
A. 
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The behavior of the fluctuating quantities is shown in Figure 3 where the 
turbulent velocity u rms and the normalized turbulent magnetic field B Tms /By 
(where -B rms is defined analog to M rms ) in dependence of the imposed magnetic 
field respective A is presented. Both quantities are significantly reduced com- 
pared to a non-magnetic but rotating convective state, indicating the trend to 
a more laminar flow. 

The ratio between the total magnetic energy and the kinetic energy is plotted 
in Fig. 3 (right). Equipartition is reached for A 0.6 indicated by the dotted 
lines. For A > 0.6 the magnetic energy clearly exceeds the kinetic energy. For 
A = 100 the magnetic energy dominates the kinetic energy by a factor of 100. 

When increasing A from 0.01 to 1000 the combined effects of rotation and 
magnetic field lead to Rossby numbers from 7 • 10~ 2 to 1 • 10~ 2 and turbulent 
Mach numbers Ma = u ims /c s from 7 • 10~ 2 to 1 • 10~ 2 . This is still much larger 
than the real Mach number in the fluid outer core which is assumed to be of 
the order of 10~ 7 . Reynolds numbers, Re = u rms d/u, are in the range from 
Re=210 for A = 0.01 and Re=32 for A = 1000. 



3.2 Patterns of flow and magnetic field 

The dynamical influence of the imposed magnetic field can be seen in figure 4 
where the z-component of the velocity near the domain faces and in a hor- 
izontal plane at z — 0.5 is shown. The left (right) panel shows a snapshot 
of the developed rotating magnetoconvection for A ps 1 (10). Upflows are vi- 
sualized in grey whereas downflows are represented in dark tones. Compared 
to non-magnetic, rotating convection the magnetoconvection pattern remains 
nearly unchanged for A<1. Many topologically not connected columnar con- 
vection cells can be seen tilted by an angle of 45° with respect to the ,2-axis 
and become aligned with the rotation axis (Taylor-Proudman theorem). 

Stronger magnetic fields lead to remarkable changes in the flow pattern. Be- 
tween A = 1 and A = 4 the quasi-regular pattern becomes more and more 
disintegrated and evolves towards a nearly two-dimensional flow as illustrated 
in the right panel of Figure 4. The convection cells are clearly elongated along 
the imposed magnetic field direction (^-direction) and show little variations of 
the convective velocity along the field lines. The sheetlike convection cells are 
again tilted inside the box and aligned with the rotation axis i? as it is the 
case for the weak-field calculations. Compared to the cases of rotating convec- 
tion or weak-field rotating magnetoconvection, the strong-field case is further 
characterized by a significant reduction of the number of convective cells. The 
nearly two- dimensionality of the flow can be explained from a condition similar 
to the Taylor-Proudman theorem for rotating spheres: For a stationary state 
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Fig. 4. Vertical velocity pattern of rotating magnetoconvection for A = 1 (left) and 
A = 10 (right). 

with small deviations from the basic unperturbed state and neglecting diffu- 
sive terms it follows from the induction equation that BidiUj = i.e. motions 
cann ot vary in the direction of the imposed magnetic field ( C handrasekhar . 



1961). 



3.3 Mean fields and a -coefficients 



Figure 5 shows the components of the horizontally averaged magnetic field in 
units of the initial field Bq for the cases A = 0.1, 1, 10. The light gray lines 
represent different snapshots within the averaging period indicating substan- 
tial fluctuations, and the thick dashed line gives the time average of these 
individual curves. There exists in all cases a significant magnetic field com- 
ponent (B x ) close to the boundaries and in the bulk having opposite signs 
in that region s. This behavior strongly differs from the results reported by 
Ziee £3 J2002) for highly stratified convection where the x-component of the 



magnetic field was found negligible ((B x ) <C (B y )). This is probab l y due to a 
much lower Taylor number employed in the simulations of IZieg le] {2OO2) and 



due to the initial two-layer configuration consisting of a convectively instable 
layer on top of a stable layer. The generation of (B x ) at the vertical bound- 
aries is a result of the perfect conductor boundary conditions for the magnetic 
field which means that no magnetic flux can cross the boundaries favoring 
concentration of magnetic flux at the top and at the bottom of the domain. 

As expected, no significant mean magnetic field in ^-direction establishes. For 
weak imposed fields local dynamo action leads to a slight amplification of (B y ) 
close to the boundaries. This effect does not take place for stronger imposed 
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Fig. 5. z-dependence of mean magnetic fields for A = 0.1, 1, 10 (from left to right). 
From top to bottom: (B x ) , (B y ) and (B z ) 

fields where, in contrast, a reduction of (B y ) occurs at the vertical boundaries. 
Usually, the dynamo a-coefficients are computed in a simplified way by 

£ x Q-xy \By) , £y &yy (-^3/) ) &z ®-zy \^y) ■ (H) 



This may not be justified with the mean field configuration obtained here, 
since the presence of a (B x ) and the presence of gradients of the mean field 
components close to the vertical boundaries lead to a contribution from non- 
diagonal coefficients of the a-tensor (describing the anisotropy) and from 
turbulent diffusivity. Their contribution cannot be easily separated without 
contradiction, so it seems more reasonable to regard the calculated EMF- 
co mponents as the mo r e rele vant quantities. This is confirmed by calculations 
of iBrandenburg et al. (1990) who showed that the expressions (11) together 
with the mean field induction equation (1) are too simple for reproducing 
the mean-field components as computed numerically from solving the full set 
of MHD-equations. However, these workers have assumed that this deviation 
occurs for shorter timescales whereas for longer timescales the correct time de- 
pendence of the mean-field is retained if a suitable quenching function for a is 
used. Due to the well-known theoretical background of characteristic proper- 
ties for the a-coefficients and to provide input data for a mean-field a 2 -dynamo 
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model we nevertheless use relationship (11) keeping in mind that it stands for 
a more or less rough approximation. 



In principle, with the present simulations it is possible to calculate the three 
coefficients a xy , a yy and a zy . Here, we want to concentrate on the first two 
coefficients a 7 



<:ri/ 



and a yy . a xy describes turbulent (radial) advection of mag- 



netic field (pumping), whereas a yy describes the production of magnetic field 
perpendicular to (B y ) via the a-effect and is therefore of profound interest. 
The vertical profile of the calculated a-coefficients is plotted in Figure 6 for 
A = 0.1, 1, 10, 100. The top row shows a xy and the bottom row shows a 



Again the grey lines represent snapshots of the coefficients at different times 
and the dashed line in each plot represents the time-averaged z-profile. 
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Fig. 6. z-dependence of a xy (top) and a yy (bottom) for A = 0.1, 1, 10, 100 (from left 
to right) at 6 = 45°. Note the reduction in scale by a factor of 5 for the A = 100 
case. 

We note that oc yy {z) shows a quite asymmetric behavior with respect to z. It is 
always negative in the lower part of the layer and positive in the upper part of 
the layer. The transition between negative and positive a-effect occurs roughly 
in the middle of the box. The peaks of Gt yy are remarkably broad leading to 
two extended zones with positive and negative a-effect having nearly equal 
amplitude. Comparable solutio ns propo r tional to sin(27rz) from quasi-linear 
calculations have been found bv ISoward (jl979h . For increasing magnetic field 
time fluctuations in the a-coefficients are obviously reduced which comes from 
the fact that the flow tends to become laminar for large A. The resulting 
vertical profiles of a yy differ from computations with stronger stratification 
where the ^-dependence is much more asymmetric. In the latter scenario the 
zero line is crossed in the lower half-box and the peak of the a-coefficients is 
higher in the upper half-box. This leads to a non-zero a-effec t when averaged 
over the entire box volume (|Riidiger and Hollerbachl . 12004). In contrast to 
this, the volume averaged a-effect is rather small in our model. 



a 



shows a drastic change in its behavior between A = 0.1 and A = 1. 
For A = 0.1 a xy has a broad peak located in the lower part of the layer and 
changes its sign close to the upper vertical boundary. For A > 1 a xy shows a 



12 



reversed z-dependence. This feature is correlated with the change of behavior 
of (By) at the boundaries between A = 0.1 and A = 1. 

Computations with a box located at 9 = 135° reveal a xy symmetric with 
respect to the equator but a yy antisymmetric (see Figure 7). The resulting 
(B) (z) closely resemble the symmetry properties of the a-coefflcients ((B x ) 
antisymmetric, (B y ) symmetric). From that it follows that the coefficient a xy 
plays the role of an (negative) advection velocity (m csc = — a xy ). The ^-profile 
of a xy shows no preferred symmetry with respect to the middle of the unstable 
layer, and a xy is - independently of A - predominantly positive. The escape 
velocity in the vertical direction is thus directed downwards for strong magnetic 
fie lds. Exactly this behavior has be en obtained in a quasi-linear approximation 
by Kitchatinov and Riidiger ( 1992f ) characterized as 'tur bulent buoyancy' . The 
same r esult has been found in nu meri cal simulations bv iDorch and Nordlund 
f)200lMOssendriiver et all (|2002h and IZieder and Riidigerl (|2003h . It can now 
be considered as a well-established phenomenon that the magnetic-induced 
turbulent pumping transports mean magnetic field downwards rather than 
upwards to the surface. 



A 



A = 10 





0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



- 1 n r 

0.0 0.2 0.4 



D.D 0.2 0.4 D.6 0.8 



V -0. 



0.0 0.2 0.4 D.6 0.8 



0.0 0.2 0.4 0.6 0.8 



D.D 0.2 0.4 0.6 0.8 



Fig. 7. z-dependence of a xy and a yy 
(B y ) (lower row) for A = 1 and A = 



(upper row) and the corresponding (B x ) and 
10 at 6 = 135° 



3.4 Kinetic helicity 



Figure 8 shows the z-profile of the kinetic helicity H^ in = ( u ' • V x it'). The 
kinetic helicity is negative in the upper part of the box and positive in the 
lower part. Comparing figure 8 with the corresponding a yy -profiles from fig- 
ure 6 the well-known relation between the signs of the a-coefficient and the 
kinetic helicity is confirmed, i.e. a ~ — l/_3_r cor (it' • V x it') with a correlation 
*mT-J~„kmm.mi Mda ESa A relation was found by 



Ossendriiver et al.l ([20011 ). Increasing the strength of the imposed magnetic 



field the change in amplitude of the helicity is roughly in accordance with the 
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Fig. 8. z-dependence of the kinetic helicity for A = 0.1, 1, 10, 100 



change in amplitude of a yy . The correlation time r cor roughly is of the order 
of 10% of the turnover time r adv 



3.5 a- quenching 



In order to estimate the quenching behavior of the a-effect we investigate 
the variation of the local maximum (minimum) value of the time averaged 
^-profile ct yy {z) in dependence of the quantity (By/B cq ) 2 where By denotes 
the initially imposed magnetic field strength and B eq is the so called equipar- 
tition field which is defined by B 2 q = fiopu 2 ms (B — >0), i.e. w rms is related to the 
non-magnetic case. In a rough approximation we have A m 6.24 • (By/ B eq ) 2 . 
The quenching curves are plotted in figure 9 where the solid (dotted) line cor- 
responds to the maximum (minimum) of a yy (z). The dashed line corresponds 
to a simple analytic quenching function of the form 

a(B) = a(B^O) -i ^. (12) 

i + 0°/^) 



A detailed analy sis of the a-quenching phenomen on for slow rotation, however, 
can be found in iRiidiger and Kitchatinov The corresponding expres- 



sions for a yy is given by the dashed-dotted line in figure 9. It overestimates the 
suppression of ct yy so that a-quenching is probably better described by equa- 
tion (12) for our case of a fast rotator. Note that the curves associated with the 
maximum respective minimum of a yy are almost identical for (By/ B eq ) 2 >0.2 
and the deviation is at most about 10% for (B®/ B eq ) 2 <0.2. 

For the dynamo number, C a = ad/rj, we obtain a value of the order of 10 
for imposed fields below (By/B cq ) 2 »s 1. This value might be large enough to 
allow global dynamo action in the sense of a 2 -dynamos but this, of course, is 
a very crude estimation. 
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Fig. 9. Quenching behavior of a yy . A value of (By/B eq ) 2 = 1 (equipartition value) 
corresponds to A ~ 6.2. The solid (dotted) line corresponds to the maximum 
(minimum) of a yy (z), the dashed line corresponds to the analytic quenching func- 
tion (12) and the dashed - dotted line represents the analytic expression derived by 
iRudiger and Kitchatinov (|l993h . 



4 Conclusions 



The main result of our computations is that the y-component of the EMF 
describing the a-effect in the azimuthal direction is positive in the upper 
part of the unstable layer and negative in the lower part on the northern 
hemisphere. We obtained the ^-profiles of the a-coefficient for a wide range of 
imposed magnetic field strength with the amplitude of the a-effect probably 
sufficient to ensure dynamo action in mean-field calculations. 

The a-effect is quenched under the influence of strong external magnetic fields 
but the suppression is not catastrophic, and significant quenching only sets in 
for magnetic fields above the equipartition value B eq . 



Some weak points in the above model calculations should be mentioned. First, 
we neglected the radial dependence of the gravitational force and the effects of 
spherical geometry, thus, the influence of curvature effects remains unknown. 
Compositional convection is ignored and the adopted parameter values are still 
far from the realistic values for the Earth. Somewhat unphysical boundary con- 
ditions have been used by assuming the conditions of a perfect conductor for 
both vertical boundaries so it cannot be ruled out that the obtained profiles 
are profoundly affected by the vertical boundary conditions. Test calculations 
with rigid boundary conditions for the velocity and / or slightly different bound- 
ar y condit i ons fo r the magnetic field as used e.g. bv lOssendriiver et al. ( 2002h 
Ziegler (2002) show, however, no major differences in the EMF profiles. It 



or 



remains, therefore, a much-promising ansatz for future mean-field calculations 
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of the geodynamo taking into account the special radial profile of the a-effect 
as it has been obtained in the present paper. 
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